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Qlh| For the lithium dimer we calculate cross sections for absorption of radi- 

ation from the vibrational-rotational levels of the ground X 1 S^" electronic 
state to the vibrational levels and continua of the excited A X E^ and B 1 !!^ 

^ | electronic states. Theoretical and experimental data are used to character- 

ed ' 
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from photoassociation spectroscopy and ultra-cold atom collision studies. The 

43 ; 

quantum-mechanical calculations are carried out for temperatures in the range 
from 1000 to 2000 K and are compared with previous calculations and mea- 
surements. 

i PACS numbers: 33.20.-t, 34.20.Mq, 52.25.Rv 
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ize the molecular properties taking advantage of knowledge recently obtained 
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I. INTRODUCTION 



The absorption spectra of pure alkali-metal vapors at temperatures of the order 1000 K 
can be a rich source of information on molecular potentials and properties. Achieving a high 
vapor pressure of lithium in experiments requires higher temperatures than the other alkali- 
metal atoms, but there are some data from heat pipe ovens |]-f§ and from a specialized 
apparatus ||]. In addition to the atomic lines the spectra exhibit gross molecular features 
attributable to transitions between bound levels of the ground electronic state and levels of 
the excited singlet states and weaker features arising from analogous triplet transitions. 

Theoretically, the envelope of the alkali-metal molecular absorption spectra can be quan- 
titatively reproduced using semi-classical models ]5|-f| and ro-vibrational structure |7j and 
continua [K| can be reproduced from quantum-mechanical models. In this paper we cal- 
culate quantum-mechanically absorption spectra for the X 1 S^-A 1 Ej and X^^-B 1 !^ 
transitions in L12. Although both semi-classical @ and quantum-mechanical |7l,|TTH calcula- 
tions have been performed and compared |7j previously for these transitions of Li 2 , recent 
improvements in the molecular data prompt the present comprehensive study. 

From photoassociation spectroscopy and cold collision studies performed in the last few 
years as well as recent theoretical work there have been significant critical tests of and 
improvements to the molecular potentials [0^0' particularly at long-range [fl5| -|l7|, and 
transition dipole moment data |~5| available for Li 2 , as well as to the value of the lifetime 
of the Li 2p state [|T~8| T p~9f| . This paper presents calculations of the spectra over the full range 
of wavelengths where absorption in the X 1 S+-A 1 S+ and X 1 E+-B 1 I1 U bands is possible. 
We calculate the satellite feature profiles at various temperatures, identify and explore the 
influence of quasibound states and the contributions of bound-bound versus bound-free 
transitions, calculate partition functions, and calculate lifetimes for the A X E+ and B 1 U U 
ro-vibrational levels. 



II. QUANTUM THEORY OF ABSORPTION CROSS SECTION 

In the quantum-mechanical formulation an absorption cross section from a vibration- 
rotation state of the lower electronic state (v", J", A") to the vibration-rotation state of the 
upper electronic state (v f , J', A') is 

rtfA-iy) = ^\(<f> M "\D(R) \<j> M ,)\ 2 g(v - ^)^ffj (1) 
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where g(y — i/y) is a line-shape function of dimension z/ -1 , Sj,,\ii is the Honl-London factor 
and z/jj = l-Bu'j'A' — £v'j"A"| is the transition frequency |2(J,^l[]. In this study, g[y — Uij) is 
approximated by 1/Az/, with Au the bin size. For a bound-free transition, the absorption 
cross section from a bound level of the lower electronic state (v"J"A") to a continuum level 
of the upper electronic state (e'J'A') can be written as 

oifAniy) = ^ K<W'A"| D(R) \^ JIA/ )\ 2 ^ff_ (2) 

where the continuum wave function 4> e ij'M is energy normalized. Free-bound or free-free 
transitions are not considered since the temperatures studied here are not high enough for 
these types of transitions to be important within the singlet manifold. 

The radial wave function can be obtained from the Schrodinger equation for the relative 
motion of the nuclei 

where V(R) is the rotationless potential energy for the electronic state, fi = 6 394.7 is the 
reduced mass of the 7 Li atoms, and E is for bound states the eigenvalue E v j\ measured 
with respect to the dissociation limit associated with the wave function <p{R) = <p v j\(R). 
Similarly, for continuum states E is the relative kinetic energy of the colliding atoms E e j\ 
associated with the energy-normalized wave function <p{R) = <f) t j\(R). 

The total absorption cross section at frequency v can be obtained by averaging over 
initial vibration-rotation levels (v",J",A") with a relevant weighting factor and summing 
over all possible bound-bound and bound-free transitions with frequencies between v and 
v + Av yielding 
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£ al^uj^J" + 1) exp[-(D e + E v „j„ K ,,)/kT\ 

U'J"v'v" 

+ E / de'ai.fj^ujn (2 J" + 1) exp[-{D e + E v „j„ K „)/kT] 

J'J"v" J 



(4) 



where ujj» is a statistical factor due to nuclear spin with the values [// (21 + 1)] = | for even 
J and [(/ + 1)/(2J + 1)] = | for odd J, for 7 Li 2 with / = |. With the zero of energy taken 
to be the potential minimum, the partition function Z\ of the lower state with dissociation 
energy D e is 

Z l = Uj»{2J" + 1) exp[-( J D e + E v „j,, K ,,)/kT], (5) 

J"v" 

assuming thermodynamic equilibrium. The resulting cross sections can be used to model 
the dimer absorption spectra in the quasistatic limit. 
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III. MOLECULAR DATA 



The adopted molecular potentials of the X 1 S^, A X S+ and B 1 !!^ states are shown in 
Fig. [TJ. The ground X X E+ state potential was constructed using the approach of Ref. [f22] . 
We adopted the recommended |23| potential obtained by Barakat et al. PH , who applied 



the Rydberg-Klein-Rees (RKR) method to measured energies. The data were connected 
smoothly to the long-range form 



V(R) = 

where V exc (R) is the exchange energy 



R* RS RIO + V "«W> 



(6) 



and the coefficients C 6 , C 8 , and Cm have been 
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calculated in Refs. [p^jp5| -[27|] and we use atomic units in this section. We adopted the 
coefficients C 6 = 1393.39, C 8 = 83 425.8 and C 10 = 7.372 x 10 6 from Ref. |§. The 
two regions were connected at R — 23.885 a® yielding a value of 8516.95 cm" 1 for the 
dissociation energy of the X 1 !!^ state of 7 Li 2 , in satisfactory agreement with the accepted 
value of 8516.61 cm" 1 | 23| . The form aexp(— bR) was used to extrapolate the potential 
at short-range where the constants were determined to smoothly connect to the innermost 
RKR points. The resulting X 1 Sj" potential yields an s-wave scattering length of 33.6 ao, 



in excellent agreement with the accepted 
assembled data. 



value of 33 ± 2 and a sensitive test of the 



For the excited A 1 ^^ state the RKR potential of Ref. [17j was adopted and smoothly 
connected at about R = 140 ao to the long-range form 



R 3 



V(R) 
-11.000 226 and C 6 



R 6 



Cs 



2 075.05 from Ref. pi and C s 



(7) 
2.705 x 



with coefficients C3 = 

10 5 from Ref. [15||. For R < 3.78 ao the data were connected to the short range form 
aexp(— bR). The dissociation energy for the A X E^ state is determined to be 9 352.194 cm -1 
in our calculation, in good agreement with the experimental value [29] of 9 352.5(6) cm -1 
given in Ref. p0| . 

The potential for the B 1 Il n state has a hump with maximum at R ~ 11 ao and various 



determinations of the hump location and height have been summarized in Refs. [31| and [30 



We adopted the IPA (Inverted Perturbation Approach) potentials from Hessel and Vidal |32 



for R < 9.35 a and the ab initio potential of Schmidt-Mink et al. |30| for 10.5 < R < 30 a 
with one additional point at R = 11.2 ao fixing the barrier maximum energy at 512 cm -1 
above dissociation. At R = 30 ao the data were connected to the long-range form of Eq. (|7|) 
by shifting down the data by 0.3 cm -1 from 10.5 < R < 30. The values of the coefficients 
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used in Eq. (g) were C 3 = 5.500 113 and C 6 = 1406.08 from Ref. |§] and C 8 = 4.756 x 10 4 



from Ref. The potential energy data for R < 9.35 a were fixed using the B 1 !^ state 
dissociation energy of 2 984.8 cm -1 , which we determined using the experimental value for 
T e of 20 436.32 cm -1 [32], the atomic asymptotic energy of 14 904.0 cm -1 , and the X X E+ 
dissociation energy of 8 516.61 cm -1 [p3fl . Finally, the data in the range 9.35 < R < 10.5 were 
smoothly connected using cubic splines. For R < 4.254 61 a the data were extrapolated 
using the short range form aexp(— bR). 

Transition dipole moments for the X 1 S+-A 1 S+ and X^^-B 1 !^ transitions are avail- 
able from ab initio calculations 0,0,0 and for X-A transitions from measured lifetimes of 
A state levels p5| . For the electronic transition dipole moments D(R), we adopted ab initio 
calculations of Ratcliff, Fish, and Konowalow |E| connected at R = 35 a to the long-range 
asymptotic form 



D^R) =D + —. 



(8) 



The value of the coefficient Dn was 3.3175 for both X-A and X-B transitions and the 



coefficient b was 283.07 or —141.53 [To] for, respectively, X-A or X-B transitions. For both 
transitions, we multiplied D^R) calculated using the above coefficients by a constant such 



that the value 0^(35) was identical to the corresponding ab initio value from Ref. |33 



to provide a smooth connection between short and long-range forms for D(R). The X-A 
dipole moment function that we adopted is consistent with that derived by Baumgartner et 
al. [^] from experimental measurements. There is no experimentally-derived dipole moment 
function for the X-B transition. 



IV. RESULTS 



Bound and continuum wave functions were calculated using the Numerov method to 
integrate Eq. (|3|). For the X and A states, eigenvalues were generally in good agreement 
with the Rydberg-Klein-Rees values used as input to the potentials constructed. [There 
is an apparent misprint for the energy of the v" = 9 level of the X 1 S+ state in Ref. |24| . 



We used 3 098.6412 cm -1 , consistent with Ref. [29"]-1 The constructed B 1 ]!^ state potential 
can reproduce the rotationless IPA energies tabulated by Hessel and Vidal typically to 
about 0.1 cm -1 , with the greatest discrepancy 0.15 cm -1 for the v' = 13 value. Calculated 
frequencies of B 1 I1 U -X 1 S+ transitions were also compared with calculations of Verma, Koch, 
and Stwalley []36| and the agreement was good, within 0.1 cm -1 . Four quasi-bound levels of 
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the B 1 II n state were found. For the rotationless potential their calculated eigenvalues are 
143.91, 276.7, 391.74 and 483.88 cm" 1 above the dissociation limit for v' = 14 to v' = 17. 

The calculated term energies of vibration-rotation states in the X X E+ state were used 
to compute partition functions using Eq. (|5]). The maximum vibrational and rotational 
quantum numbers in our calculations are 41 and 123, respectively, for the X 1 S+ state. 
In the harmonic approximation for ro-vibrational energies, the partition function can be 
calculated using the simple expression 

123 

Zi = Z R Z V « ]T ujjn(2J" + 1) exp[-hcB e J"(J" + l)/kT] 

J"=0 
41 

x J2 exp[-hu e {v" + l/2)/kT], (9) 

v"=0 

which assumes that the term energy can be described by the first terms of the power se- 
ries with respect to vibrational and rotational quantum numbers. Using constants v e jc = 
351.3904 cm -1 and B e = 0.672 566 cm" 1 the partition function from Eq. (^j) was calcu- 
lated and it is compared with the partition function calculated from Eq. @ for the X X E+ 
state as a function of temperature in Fig. |2|. The anharmonicity of the potential for higher 
vibrational levels accounts for the differences between the two results with increasing tem- 
perature. For J > 2 the X 1 E+ state supports quasibound vibrational levels. The expression 
Eq. dB|) for Zi does not specify whether quasibound states are to be included or not in the 
summations. We evaluated Z[ with and without the quasibound levels to ascertain their 
importance and the results are shown in Fig. |2|. The effect of the additional levels becomes 
increasingly significant with higher temperature. For the present study covering tempera- 
tures between 1000 and 2000 K there is not a significant distinction between Zj, Zi, and the 
result with the inclusion of the quasibound states. 

The molecular fraction can be calculated using the expression 

[N Lh ]/[N u } 2 = (QujQh) eMDe/kT), (10) 

where the atomic partition function Qu is 2(2nmukT / h 2 ) 3 ^ 2 , with the electronic partition 
function for the atom well-approximated by the spin degeneracy of 2 for the temperatures 
studied in the present work, and the molecular partition function Qu 2 is (2nmu 2 kT / h 2 ) 3 ^ 2 Zi, 
with the electronic partition function for the X X E+ state taking the value 1. The molecular 
fraction Eq. (|IUD is plotted in Fig. [| The absorption coefficient k(u) can be obtained if the 
atomic density is known from 

k{y) = [NuM")- (11) 
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A. Lifetimes 



Lifetimes of the various ro-vibrational levels of the A x £+ state were measured JJS] and 



calculated, see for example, Refs. [IT^.BRBl] • We calculated spontaneous emission transition 



probabilities and lifetimes of rotational- vibrational levels of the A x £+ state in order to test 
the adopted transition dipole moment. The spontaneous emission rate from a bound state 
(v'J'A') to a bound state (v"J'A") is 

A{v'J'A'- v"J'h") = -^g \ ((j)v"j'A" \ D(R) l^j^f , (12) 
where the electronic state degeneracy is 

g= (2-W) (13) 

* — <J0,A' 

and we have neglected change in the rotational quantum number. The total spontaneous 
emission rate from the upper level (v'J'A') can be obtained by summing over all possible 
transitions to bound and continuum states 

Aiy'J'A') = Y, A(v'J'A'; v" J' A") + £ / de"A(v'J'A'; e"J'A"), (14) 

v"A" A" J 

where A(v'J'A'; e"J'A") is the spontaneous emission probability to a continuum energy e" 
with partial wave J'. The lifetime is 

r = l/A(v'J'A'). (15) 



Lifetimes of levels v'J' of the A x £+ state were measured by Baumgartner et al. ||35|| . The 



A X E+ state is affected by indirect predissociation via the a 3 £+ and 1 3 IT+ states |39],f|0 



The measured lifetimes r m of vibration-rotation levels thought to be unaffected by indirect 



predissociation taken from Ref. [35 and corresponding calculated lifetimes r c are presented 



in Fig. |3] along with calculated term energies expressed relative to the potential minimum of 



the A state. The energies are plotted in the order of the values listed in Table 1 of Ref. |35 
and correspond to a range of values of (v', J') from (0,15) to (24,25). The agreement between 
r c and r m is good within the experimental precision of ±2 percent P5fl . The quasi-bound 
levels and continua of the X X E+ state inside or above the high centrifugal potential barriers 
are found to be important in calculating lifetimes of high J levels. For instance, by including 
transitions to three quasi-bound levels and the continuum states the calculated lifetime of the 
A X E+ [y' = 20, J' = 50) level changes from 22.7 ns to 19.3 ns which is close to the measured 
lifetime of 18.66 ±0.37 ns []3"5|]. The calculated and measured lifetimes agree well up to about 
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4 500 cm 1 in agreement with the findings of Ref . |30| . From approximately 5 000 cm 1 and 



higher the experimental lifetimes slightly decrease relative to theory by about 0.8 ns or five 
percent, as demonstrated in Fig. |3|. We investigated whether the 2 1 £+ state might supply 
an additional spontaneous decay channel, but the theoretical value T e = 20 128 cm -1 [30] for 



this state appears to place its minimum at around 6 000 cm -1 relative to the minimum of the 
A X E+ state, apparently ruling this mechanism out. The reason for the significant downturn 
of experimental lifetimes for higher term energies is currently not understood; however, the 
overall excellent agreement between our calculated lifetimes and the measurements gives 
us confidence in our molecular data and calculational procedures. We compare selected 
examples from the present results with calculations by Watson |37] and by Sangfelt et al. [38 



in Table |I[ The calculations of Sangfelt et al. are larger than ours, probably because 
theoretical transition energies were used. A simple rescaling using experimental energies, as 
pointed out by those authors, yields lifetimes in good agreement with the present work. The 



dipole moment function calculated by Watson 3]J is in good agreement with that adopted 
in the present study and cannot account for the shorter lifetimes obtained in that study. We 
present a more extensive tabulation of A X S+ lifetimes in Table [TJ covering the same values 



tabulated in Table VII of Ref. 38 



Lifetimes for levels of the B 1 Ii u state have been calculated by Uzer, Watson, and Dal- 



garno |4"I |, Uzer and Dalgarno j39|], and Sangfelt et al. |38|] and there appear to be no 
experimental data. In Table |TTT| we compare our calculated lifetimes of selected B l U u lev- 
els with available calculations for higher values of J'. The present results lie between the 



calculations of Uzer et al. an d those of Sangfelt et al. [pq| . The lifetimes calculated by 
Uzer et al. are longer than ours because their transition dipole moment function, calculated 
using a model potential method, is smaller than the ab initio dipole moment of Ratcliff 
et al. adopted in the present study. The dipole moment function calculated by Sangfelt 
et al., on the other hand, is in good agreement with that adopted in the present work. As 
those authors pointed out, and as our results illustrate, the utilization in their calculation of 
calculated excitation energies which were larger than experimental energies yielded lifetimes 
that were too short. In Table [IV] we present a more extensive tabulation of B 1 II M lifetimes 
covering the same values tabulated in Table VIII of Ref. PS . 



B. Absorption cross sections 

Absorption spectra arising from molecular singlet transitions at the far wings of the 
atomic 2p line at 671 nm consist of a blue wing due to X 1 E+-B 1 n u transitions and a 
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red wing due to X 1 S+-A 1 S+ transitions. Calculations for bound-bound (bb) and bound- 
free (bf) absorption cross sections were carried out separately using Eq. with Au = 
8 cm -1 . The results for the total (the sum of bb and bf) absorption cross sections at 
temperatures of 1000 K, 1500 K and 2033 K are given in Figs. f|-|6|. The ratios of the 
peak cross sections of the X 1 S+-B 1 I1 U wing to those of the X 1 S+-A 1 S+ wing are higher 
at lower temperatures. As temperature increases, absorption spectra spread out from the 
peak spectral regions and there emerges near 900 nm a satellite feature arising from the 
|42l in the X 1 S+-A 1 S+ difference potential and the maximum of the transition 



minimum 



dipole moment function We also show the bf contributions to the cross sections 

on each plot. The bf component contributes mainly to the extreme blue part of the blue 
wing and increases in magnitude significantly as the temperature increases. It is found that 
transitions to quasibound and continuum levels of the B l U u state contribute significantly 
to the total absorption spectra in the case of X 1 T,'^-'B 1 Il u transitions, apparently because 
there is less vibrational oscillator strength density in the discrete part of the spectrum for 
the B state compared to the A state. Transitions into quasibound states of the A X E+ or 
B l U u states have been included in the results for the total cross sections in Figs. f|-|6|. 

Theoretical quantum-mechanical calculations for absorption cross sections from the 
X X S+ state to the A X S+ state over the spectral range 600-950 nm were carried out by 
Lam et al. [[5J using a constant dipole moment of 6.58 D at 1020 K, for which bf transitions 
are not significant. Our calculations in Fig. |] are about a factor of 10 less than the result 
shown in Fig. 5 of Ref . |7| , but agree well both in overall shape and in details of finer struc- 
tures. We repeated the calculations using the constant dipole moment of Ref. for both 
classical and quantum-mechanical cross sections and although these two results agreed with 
each other, they were also a factor of 10 less than the result shown in Fig. 5 of Ref. J7j. 
Thus it appears to us that there may be a mislabeling of the vertical axis in Fig. 5 of Lam et 
al. [A similar calculation that we performed {yj for Na2 at 800 K is in complete agreement 
with Fig. 4 of Ref. 0.] 

Calculations of absorption spectra at 2033 K over the spectral range 450-750 nm pre- 
sented in Fig. ^| are in good agreement with the measured values of Erdman et al. Q and 
with quantum-mechanical calculations performed by Mills [JCTf . The experimental study of 
Erdman et al. [|J] involved an investigation of molecular triplet states [Q and did not explore 
the satellite feature at 900 nm. The calculations over the range 450-750 nm by Mills [TT| 



included triplet molecular transitions and are not directly comparable with the present re- 
sults. Nevertheless since the singlet transitions dominate the absorption we find excellent 
qualitative agreement with the calculations presented by Mills. 
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TABLES 

TABLE I. Comparison of calculated lifetimes in ns for ro-vibrational levels of the A 1 ^^ state. 



J' 



Watson [37 1 



Sangfeldt et al. [38| 



This work 



7 
9 

20 



15 
5 



16.8 
16.9 
17.3 



19.25 
19.29 



18.55 
18.65 
19.04 



TABLE II. Lifetimes in ns for ro-vibrational levels of the A X S^ state calculated as described 
in the text. 



v' 


J' = 


J' = 9 


J' = 15 





17.74 


17.77 


17.82 


1 


17.87 


17.90 


17.94 


2 


17.98 


18.01 


18.06 


3 


18.09 


18.12 


18.17 


4 


18.20 


18.23 


18.27 


5 


18.30 


18.33 


18.37 


6 


18.39 


18.42 


18.46 


7 


18.48 


18.51 


18.55 
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TABLE III. Comparison of calculated lifetimes in ns for ro-vibrational levels of the B 1 U U state. 





5 



J' 
15 
9 



Uzer et al. (41 ] 
8.3 
8.5 



Sangfeldt et al. @ 



6.83 
7.20 



This work 
7.66 
7.95 



TABLE IV. Lifetimes in ns for vibrational-rotational levels of the B 1 U U state calculated as 



described in the text. 



v' 


J' = 1 


J' = 9 


J' = 15 





7.65 


7.65 


7.66 


1 


7.70 


7.71 


7.72 


2 


7.76 


7.76 


7.77 


3 


7.81 


7.82 


7.83 


4 


7.88 


7.89 


7.90 


5 


7.94 


7.95 


7.97 


6 


8.02 


8.03 


8.04 


7 


8.10 


8.10 


8.12 
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FIGURES 




5 10 15 20 25 
internuclear distance (a ) 

FIG. 1. Adopted potentials for the X 1 !^, A X S+, and B 1 ]!^ electronic states. 
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temperature (K) 

FIG. 2. Comparison at various temperatures of the partition functions Zj, calculated using 
Eq. (|9|) and experimentally determined spectroscopic constants, curve A, and Zi, from Eq. (|5|) 
and numerically determined eigenvalues, curves B and C. Inclusion of quasibound states in the 
calculation of Zi results in curve C as discussed in the text. Curve D represents the molecular 
fraction [-ZVLi 2 ]/[-/VLi] 2 5 Eq. (|i~0|), as a function of temperature. 
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FIG. 3. Comparison of our calculated lifetimes (diamonds) and measured lifetimes (circles) 
from Ref. [pSfl . The error bars indicate the quoted experimental uncertainty of ±2 percent. The 
levels given are those that were measured, ordered by increasing energy, as listed in Table 1 of 
Ref. fl35l. 
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FIG. 4. Total absorption cross sections from X 1 S+~A 1 S+ and X^^-B 1 ]^ transitions in- 
cluding bound to bound and bound to free transitions at a temperature of 1000 K. The satellite 
feature near 900 nm does not appear at this temperature and bound-free absorption (dotted curve) 
is insignificant. The inset presents a magnified view of the bound-free contribution. 
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400 500 600 700 800 900 1000 

wavelength (nm) 

FIG. 5. Total absorption cross sections from X 1 S+-A 1 S+ and X 1 S+-B 1 II U transitions 
including bound to bound and bound to free transitions at a temperature of 1500 K. As the 
temperature increases, the absorption spectra are distributed over a wider spectral range and the 
ratio of the peak cross sections between X 1 T,+~B 1 U U and X^^-A 1 !!^ bands decreases. The 
satellite feature near 900 nm and bound-free absorption (dotted curve) are noticeable at this 
temperature. The inset presents a magnified view of the bound-free contribution. 
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FIG. 6. Total absorption cross sections from X 1 E+-A 1 E+ and X 1 E+-B 1 II U transitions 
including bound to bound and bound to free transitions at a temperature of 2033 K. At this 
temperature the satellite feature near 900 nm is now prominent and bound-free absorption (dotted 
curve) contributes significantly at 450 nm. The inset presents a magnified view of the bound-free 
contribution. 
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